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In order to discuss the effectiveness of inverter performance, it is essential to 
concentrate on its losses. Thus, knowing the junction temperature (TJ) value 
of the inverter switches insulated-gate bipolar transistor (IGBTs) under 
different loading conditions allows determining their reliability and 
designing the right heat sink. The old methods to calculate the TJ are 
inaccurate, difficult to implement, and require mathematical analysis. So, in 
this work, a new approach was used to calculate the TJ in a simplified way, 
using the SIMSCAPE environment in the MATLAB simulation. The IGBT 
in the SIMSCAPE environment is of an improved design and has a direct 
thermal port. The TJ estimation is more accurate and it relies only on the 
thermal impedance model, so we do not need to compute the power losses as 
was done in previous studies. This functionality is only accessible in more 
recent versions of MATLAB (2019-2023). With different motor load 
torques, the IGBT junction temperature is evaluated, and the effect of the 


selected heat sink on the value of TJ is evaluated when it is under the effect 
of natural cooling and air-forced cooling. 
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1. INTRODUCTION 

Power electronics (PE) devices are widely used in industrial and household settings in contemporary 
civilization for purposes of controlling and producing electricity [1]. Power transistors, particularly insulated- 
gate bipolar transistor (IGBT) modules, are in control of the power converter industry thanks to their superior 
performance, inexpensive cost, dependability, lightweight, wind power generation, and motor drives [2]. The 
semiconductor devices' operating junction temperature has a significant influence on IGBT power modules' 
lifespan [3]. One of the most important factors to take into account is the junction temperature, which is 
directly related to inverter power loss and overheating. These elements can lead to wear and tear as well as 
failures. The semiconductor chip junction temperature prediction is necessary to achieve accurate 
measurement and avoid chip burning [4]. Many approaches to simulating the junction temperature were 
previously in place, such as finite element simulation and three-dimensional modeling. while the two 
approaches above make it impossible to extract the internal and structural parameters [5], [6]. In addition to 
the most common method, which is to make a two-part simulation through coupling, the RC network model 
and loss calculation model, which are built on the foundation of a data sheet, are then electro-thermally 
coupled, with the relevant parameters of the inverter system being entered into the loss calculation model to 
obtain the thermoelectric simulation [7]. In this study, the junction temperatures of all IGBTs can be directly 
estimated by using the SIMSCAPE environment in MATLAB/Simulink and building an impedance thermal 
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model alone based on the manufacturer's data. This work focused on overcoming these restrictions and 
developing a junction temperature estimating system that included various load patterns for induction motors 
and types of applied coolant systems. 


2. METHOD 
2.1. Thermal resistance model 

As inverters have been developed to drive induction motors and supply important applications, it 
has become crucial to understand how to regulate the IGBT junction temperature, which is the power module 
of an inverter. The IGBTs are often restricted by their junction temperature, and if the junction temperature is 
not correctly controlled, there is a significant risk of failure due to overheating [8]. The IGBTs are 
responsible for all power loss; the heat produced by the IGBT is mostly transferred to the substrate, heat sink, 
and air through heat conduction. As a result, the thermal resistance of the IGBT module's entire heat transfer 
process can be divided into only three components. The IGBT chip's thermal resistance to the heat sink is 
initially measured as its contact thermal resistance, followed by the heat sink's solid thermal resistance and 
the heat sink's convective heat transfer resistance to the ambient [9], [10]. 

Actually, the heat dissipation route diagram's shows that the thermal resistance between the chip and 
the heatsink is the total thermal resistance of each layer in series [11], as illustrated in (1). 


RJA = RJC + RCH + RHA (1) 
The following equation represents each layer's thermal resistance and heat capacity: 


Rth =— (2) 


L 
kA 
Cth = cpAL (3) 


where, Rth is the thermal resistance per each layer, L is its thickness, k is its thermal conductivity, A is its 
area for dissipating heat, Cth is the material's total heat capacity, c is its specific heat capacity, and r is its 
density. For an accurate TJ estimate and a consistent conformance of its junction temperature, the IGBT 
module in the three-phase inverter employs the junction temperature estimation approach. This task can be 
achieved by simulating the thermal network model (Foster model) using MATLAB/Simulink SIMSCAPE 
environment, which were created with the manufacturers document in [12]. Figure 1 represent the IGBT and 
the freewheeling diode module's complete heat transfer process in the thermal resistor network [13]. 
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Figure 1. Heat transfer mechanism 


2.2. Heat sink characteristics 

To keep the device junction temperature at a low-level during power flow through it, it is important 
to select a heat sink in accordance with chip losses [14]. Conduction and convection-type heat transfer are 
used by the majority of electronic device cooling systems for thermal management [15]. The heat sink 
material and size are taken into consideration when the amount of power dissipation in the device rises, and 
the material's thermal conductivity should be as high as feasible. The finned extruded aluminum form is the 
most often utilized kind of heat sink that is commercially available, as shown in Figure 2. 

This fundamental radiator form is best for forced air and natural convection cooling because a large 
surface area is exposed to the surrounding air inside a small area. The capacity of aluminum to transfer heat 
to the ambient is measured in °C/W of power dissipation to choose the appropriate kind of heat sink, consider 
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the power dissipation and thermal data from the datasheet provided by the device manufacturer, to select the 
heat sink, the following parameters are required: target junction temperature (TJ), maximum ambient 
temperature (Ta), thermal resistance value Rth(ja), and heat sink thermal resistance value Rth(sa) [16]. 
Power chips like IGBT often employ forced air cooling with a heat sink and fans to dissipate heat or liquid 
cooling techniques because natural convection can no longer dissipate heat due to the fast growth of IGBT 
modules [17]. Figure 3 demonstrates the change in the thermal resistance value of the forced-air heat sink [18]. 
The value of Rth(ja) and value of Rth(sa) can be obtained from (4) and (5). 


Rth(ja) = (Tj — Ta)/Pdiss (4) 


Rth(sa) = Rth(ja) — Rth(jc) — Rth(cs) 
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Figure 2. The common finned heat sink Figure 3. Heat dissipation of an IGBT 


2.3. Three-phase voltage source inverter fed induction motor 

Three-phase induction motors (IM) are often the most used in the industry since need less 
maintenance and operate with great reliability [19]. Moreover, customers are increasingly interested in highly 
efficient electric motors since the price of electric energy has increased in recent years [20]. Figure 4 shows a 
capacitor-DC link for a typical three-phase voltage-source inverter (VSI) with a six-switch connected to the 
three-phase IM. The three-phase voltage source inverter consists of six semiconductor switches (S1 to S6) 
that conduct in a 180° phase shift between their upper and lower switches [21]. Six operational periods are 
separated by 60 degrees in VSI. For the inverter to operate symmetrically, one upper switch, two lower 
switches, or two upper switches and one lower switch are tuned on for each interval. There are a total of eight 
different combinations in which the three-phase VSI may operate. These are shown by the switching 
positions on the top or lower side [22]. 

From the dynamic representation of the IM model, the three-phase currents can be presented in a 
three-axis coordinate. The stator current and the angular frequency, which are described in the three-phase 
coordinate frame as follows [23]: 


isa = Im sin(wo * t) (6) 
isp = Im sin (wo *t+ =) (7) 
isc = Im sin (wo *t+ =) (8) 


transforming the stator current from three-phase coordinate to two-phase coordinate as in [23]: 


is = = (ig +a * ip + a? * ic) (9) 
a= eT = (4S (10) 
a? sels = -S (11) 


in addition, to transform the other variables, such as electromagnetic torque (Te), the same transformation of 
coordinates is used [24]: 
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Vy = Rs * is +2 + jws + Y, (12) 
0 = R, +i, + + j(ws — w) «V, (13) 
V, = Ls xi; + Lm *i, (14) 
Y, = Lr x i, + Lm * i, (15) 
T, = =* p * Re{W5 * is} = — «p * Re(r x ir} (16) 


2.4. Space vector PWM technique 

Numerous industrial applications, including energy savers, motor soft starters, are utilizing different 
voltage modulation methods. Inverters controlled by support vector machine (SVM) provide the lowest 
harmonics and highest efficiency when compared to other modulation techniques [25]. Improvements in the 
quality and consistency of the drive circuit are becoming more and more necessary as a result of ongoing 
research on improving the performance of motor drives. The development of power semiconductor devices 
has greatly increased the power handling capacity of IGBT switches. SVM inverters are utilized extremely 
efficiently and affordably in low, medium, and certain high-power applications [26]. The SVPWM and its 
switching operation is presented in Figure 5 [27]. 


¥,(010) Y,(110) 
x d 


Figure 4. Three phase inverter fed IM Figure 5. Space vector PWM diagram and switching operation 


SVPWM samples the reference vector V* at a frequency fs (Ts=1/fs). The alpha-beta transformation 
is used to create the reference vector V* from three-phase 120? phase shift references. 


Ya = Van — 1 Vbn — 2 Vcn (17) 


v3 v3 


VB = Van — 7 Vbn — Vc (18) 


The magnitude of the reference voltage is given by (19). 
Vref = J (Va)? + (VB )? (19) 


There are six distinct sectors spaced 60 degrees apart on the af -plane, the unitary ramp is generated using 
the ramp generator at the PWM switching frequency, the switching sequence uses this ramp as its time basis, 
when determining when to apply a voltage vector to a motor, the switching time calculator is employed, the 
block input is the sector in which the voltage vector lies, the switching time calculator and the ramp generator 
provide the gate logic with the timing sequence and ramp, respectively, in order to properly activate the 
inverter switches, this block compares the ramp and gate timing signals [28]. 


3. RESULTS AND DISCUSSION 

The MATLAB/Simulink 2022 package is used with the SIMSCAPE environment to run the 
simulation. three-phase inverter with an induction motor load (IM). Figures 6 and 7 show the components of 
the circuit model used to simulate the proposed method. The system parameters of the simulation model were 
listed as follows: rated power is 7 KW, rated speed is 1800 rpm, and p=4. The switching frequency is 5 kHz, 
the fundamental frequency is 60 Hz, the DC-link voltage is 400 volts, the ambient temperature is 25 ?C, and 
the IGBT's maximum temperature is 150 °C. The cooling system's fan speed is 1000 ft/min. 
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Figure 6. Block diagram of the overall system 
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Figure 7. Thermal model of an IGBT with heat sink in SIMSCAPE environment 
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In this study, four load torque values (10, 15, 20, 25) Nm are used of this system to monitor changes 
in the junction temperature. It is found that the value of the junction temperature increases as the load torque 
rises; the maximum value of TJ was recorded at a torque of 10 Nm is (64.5 °C), it recorded (71.7 °C) at a 
torque of 15 Nm, it recorded (77.24 ?C) at a torque of 20 Nm and recorded (82.46 ?C) at a torque of 25 Nm. 
The maximum allowable TJ provided by the manufacturer is 150 °C. As illustrated in Figures 8 to 11. Figure 12 
shows how the value of the junction temperature changes under the natural cooling of the heat sink when an 
additional load torque is applied suddenly, where the value of TJ is 83 °C. Figure 13 shows how the value of 
the junction temperature changes under forced air cooling of the heat sink when an additional load torque is 
suddenly applied, where the value of TJ recorded is 61.2 °C. Figure 14 shows the phase current and line to 
line voltage at load torque of 15 Nm. 
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Figure 8. Junction temperature at 10 Nm 
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Figure 9. Junction temperature at 15 Nm 
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Figure 10. Junction temperature at 20 Nm 
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Figure 11. Junction temperature at 25 Nm 


Int J Pow Elec & Dri Syst, Vol. 14, No. 2, June 2023: 653-661 


Int J Pow Elec & Dri Syst ISSN: 2088-8694 Oo 659 


Torque 


Torque (Nm) 


Time (seconds) 
TJ, Tmax 
160 
—~140 
Q 
2120 
2 
100 
E 
È 8 
S 
S 60 
c 
2 
Er 405 
20 
0 1 2 3 4 5 6 7 8 9 10 
Time (seconds) 


Figure 12. Junction temperature under natural cooling 
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Figure 13. Junction temperature under air forced cooling 
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Figure 14. Line voltage Vab and phase current Ian 


4. CONCLUSION 

Calculating the junction temperature for IGBT is very important because it gives a predictive 
indication of the power converter's efficiency. In this work, the SIMSCAPE environment was used in the 
MATLAB/Simulink simulation to calculate IGBT junction temperature. as a new proposed physical 
modeling approach. Based on this new and direct modeling, it is possible to easily obtain the value of TJ by 
analyzing the equivalent heat-transfer model of the IGBT module. This can be achieved only based on the 
manufacturer's thermal impedance information. Analytical results point to the influence of the load torque on 
the value of the junction temperature. The results show that the TJ increased with the increase in IM load 
torque. In addition, the results show that under the same working conditions of the system, but by using 
forced air cooling of the IGBT instead of natural cooling, the value of the junction temperature is reduced by 
22 °C, and as a result. The performance of the IGBT will improve and its efficiency will increase while 
remaining within the allowable temperature limit. 
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